Screening ineffectiveness and THz emission at bare LO phonon frequencies 
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Within a hydrodynamic approach we investigate the dynamics of an inhomogeneous electron-hole 
gas coupled to phonons in Te and the corresponding THz emission. We find that the longitudinal 
inhomogeneity of the photogenerated electron-hole gas - due to short absorption lengths in Te - 
gives rise to screening ineffectiveness for non-zero wave- vector modes. This allows for THz dynamics 
and emission at the bare LO phonon frequency ulo even at high carrier densities. Screening 
ineffectiveness thus provides an appealingly simple explanation for the existence of bare modes at 
vlo in longitudinally inhomogeneous systems such as Te; no lateral inhomogeneity of the excitation 
spot is needed here. 



Terahertz-range dynamics of non-equilibrium carri- 
ers interacting with phonons- is a very active field 
of research in semiconductors!!! Pump-and-probe time- 
resolved spectroscopies with ultrashort laser pulses have 
played a crucial role in unveiling many interesting fea- 
tures in these complex systems. The appearance of 
time-dependent electric fields following the pump excita- 
tion - due to screening of depletion-layer fields, charge- 
separation fields, lattice polarization fields, etc. - gives 
rise to electro-optic induced reflectivity/transmittance 
changes and terahertz emission due to dipole radiation. 
Experimentalists probe the system in essentially two 
ways: by measuring (i) the reflectivity/transmittance 
changes of a probe pulse and/or (ii) the emitted radi- 
ation. These canonical experimental geometries yield, 
in principle, complementary information on the system 
characteristic frequencies. 

In n-doped GaAs, for instance, the usual two branches 
L + and L_ describing the long wavelength plasmon- 
phonon coupled modes as a function of the carrier den- 
sity are easily, mapped out from electro-optic modulation 
experiments!] As pointed out in Ref. ||, the L + mode is 
usually interpreted as a bare LO phonon mode v^o at low 
carrier densities, while the L_ mode is referred to as a 
screened LO phonon mode vto at high carrier densities. 
Interestingly enough, beating of these TO and LO modes 
is observed in the experiments. This is due to the simul- 
taneous excitation of these modes due to the lateral inho- 
mogeneity of the pump pulse which generates regions of 
high and low carrier densities - sampled simultaneously 
by the probe pulse. In contrast, terahertz emission spec- 
troscopy in GaAs provides information only on the L_ 
branch at low densities; L + and LO phonon frequencies 
(> 8 THz) are not accessible with the current spectral 
range (< 7 THz) of the state-of-the-art antennas.ta 

In Tellurium on the other hand, both reflectivity and 
THz emission measurements can probe coupled-mode fre- 
quencies. These experiments are performed in both the 
c || and cj_ geometry. These denote the either parallel or 
perpendicular orientation of the c axis of the hexagonal 
lattice of Te with respect to the excitation surface. Oscil- 
lations at bars-,^LO are seen in both experiments for the 
cm geometries^ In the cj_ geometry only THz emission 



is possible; a peak at v>lo is also evident.i,0 Similarly to 
GaAs, the existence of these peaks in the Te response can 
in principle be explained via the lateral inhomogeneity of 
the excitation spot. There is, however, a more appealing 
explanation we propose in this work. 

In this letter we theoretically investigate the coupled 
dynamics of electron-hole-phonon modes in Te and the 
corresponding THz emission. Our semiclassical descrip- 
tion incorporates phonons and its coupling to carriers 
via polarization charges in Poisson's equation. Our main 
finding is the observation that the intrinsic longitudi- 
nal inhomogeneity of the photogenerated electron-hole 
plasma - due to the short absorption lengths in Te - 
leads to screening ineffectiveness for non-zero wave vec- 
tors. Screening ineffectiveness allows for dynamics and 
terahertz emission at bare LO phonon frequencies - even 
at high carrier densities - thus providing a possible alter- 
nate/complementary explanation for the observed bare 
LO modes in the THz response of Te. We suggest a 
simple experiment to tell apart lateral- (q = 0) and 
longitudinal-inhomogeneity (q ^ 0) contributions. A di- 
electric function analysis (Debye Hiickel and RPA) cor- 
roborates our numerical findings. 

Ingredients. Any reasonable attempt to describe 
pump-and-probe THz emission experiments has to nec- 
essarily address the coupled-field dynamics of electrons, 
holes, and phonons. Our semiclassical approach com- 
bines (i) moments of the Boltzmann equation describ- 
ing charge and current densities, (ii) a phenomcnolog- 
ical driven-harmonic-oscillator equation simulating LO 
phonons, and (iii) Poisson's equation providing electron- 
and hole-phonon coupling via polarization charges. 

Model. Our one-dimcnsionalH model allows motion 
only perpendicular to the relevant semiconductor sur- 
face being excited. We consider carriers with isotropic 
masses m* = 0.06m and m* h = 0.114mo (and charges q e 
and qh) in parabolic bands.u We assume the electrons and 
holes are thermalized to._distmct constant temperatures 
T e and Th, respectively ;E3 these are, however, different 
from the lattice temperature. By considering the semi- 
classical Boltzmann equation for the distribution func- 
tion within the usual relaxation-time approximation, we 
derive equations for both particle and particle current 
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densities of electrons (e) and holes (h), n e (h)(x,t) and 
J e (h)(x,t) = n e{h) (x,t)v e{h) (x,t), respectively; v e{h) (x,t) 
denotes the electron (hole) mean velocity. These are 



dn e ( h) (x,t) dJ e(h) (x,t) 



dt 



Ox 



= G(x,t), 



(1) 



Emitted THz radiation. The longitudinal charge oscil- 
lations in the system are essentially electric dipoles and 
thus emit electromagnetic waves. Within the dipole ap- 
proximation, the emitted signal is proportional to the 
second time derivative of the internal electric field. Wc 
define the radiated field E rac [ by 
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where G(x,t) is the generation term (laser). The elec- 
tric field E(x,t) in Eq. (||) is calculated via Poisson's 
equation 



dE 



dx eqc 



1 , , x ^Tov /£ o(eo - £oo) dw 

{qhrih + qerie) -5-, (3) 



dx 



which accounts for both electron and hole charge densi- 
ties and lattice polarization charges ~ dw(x,t)/dx. We 
assume the vibrational coordinate w(x,t) in Poisson's 
equation obeys 



1 9W 2 / 7 Crn ri\ 

dt 2 '' ~~dt + uj to w = w TOV £ o(eo - eoo)E, (4) 

which is driven by the total electric field due to electron- 
hole charge separation and phonon polarization charges. 
Note that Eq. (Q) describes longitudinal oscillations. Wc 
account for damping of electrons, holes and phonons via 
proper relaxation times [r e , and r p h, respectively]. 

Equation (Q) is the continuity equation for electrons 
or holes; laser pumping of carriers is accounted for by 
the generation term (gaussian in time and exponen- 
tially attenuated in x). Equation (J2J) states momen- 
tum conservation [the momentum density is Ve(h) { x i t) — 
n e ( h ){x,t)m* e{h) v e ( h) {x,t)]. Implicit in (g) is the as- 
sumption of thermalized carriers with average energies 
(ux)e(h) = \k B T e{h] + ^m* e(h) v 2 {h) (x,t). With this as- 
sumption we are essentially truncating the sequence of 
equations for higher-order moments. By coupling Eqs. 
(Q) and ([|) , we extend the phenomenological description 
of Born and HuangEJ in that we allow the internal elec- 
tric field to drive lattice vibrations and vice versa. Our 
description here generalizes the previous investigations 
of Kuznetsov and Stantonu who treat homogeneous sys- 
tems and introduce averaging over the spot profile and 
Dekorsy et alU who use a drift-diffusion type model. 

Numerical solution. The set of six coupled differential 
equations discussed above is tricky to solve numerically. 
By carefully using standard explicit schemes for the prop- 
agation of partial differential equations,L3 we are able to 
obtain n e (x,t), J e (x,t), rih(x,t), Jh(x,t), E{x,t), and 
w(x,t) in a wide range of system parameters. As we dis- 
cuss next, both the radiated signal and the reflectivity 
response are related to E(x, t). 



which accounts for oscillating dipoles, weighed by the ab- 
sorption profile, from within only the relevant absorption 
length a -1 . Note that the radiated field E ra( i is p polar- 
ized. 

Figures 1(a) and 1(b) show the calculated emitted sig- 
nals and their Fourier transform, respectively, for several 
excitation densities n exci at fixed electron and hole tem- 
peratures, T e — 2000 K and Th = 1050 K [see also inset 
in Fig. 1(c)]. The time-dependent signals 1(a) show an 
initial rise and a dip followed by small oscillations [simi- 
lar to the experimental detected signal in the inset]. As 
shown in 1(b), these small oscillations have frequencies 
close to that of the LO phonon in Te, vlo ~ 2.82 THz 
(c_l geometry); the broad-band peak in 1(b) is due to 
electron-hole charge separation as discussed in |(| 

The detected theoretical signals in 1(c) are convolutions 
of the emitted signals in 1(b) and the antenna response 
[dotted line in 1(b)]. The convoluted signals 1(c) are 
very similar to the recent data in 1(d) from Ref. 0. In 
particular both theoretical and experimental spectra dis- 
play a broad-band peak followed by a dip at the screened 
LO frequency vto ^2.6 THz and a subsequent peak at 
the bare LO phonon frequency vlo ~ 2.82 THz. Note 
that the signals get stronger with increasing excitation 
densities [curves (1) - (6)]. This intensity enhancement 
of the signals follows from the amplitude increase of the 
internal fields with n exc . While the spectral features at 
the phonon frequencies hardly move with increasing exci- 
tation densities, the broad-band features slightly shift to 
higher frequencies. This shift is due to the faster electron- 
hole charge dynamics at higher densities!] The relative 
intensities of the bare LO peak and that of the broad- 
band peak strongly depend on the antenna response. 

Antenna response. We assumed a simplistic antenna 
response with a high sensitivity around 0.5 THz dropping 
essentially to zero above 3 THz [dotted line in 1(b)]. This 
response weighs up the contribution of the broad-band 
peak in contrast to that of the bare LO phonon. The de- 
tected/convoluted signals 1(c) are then "distorted" thus 
displaying weaker LO phonon peaks [cf. 1(b) and 1(c)]. 

Normal mode analysis. By linearizing our equations 
and neglecting damping for simplicity, we obtain the sys- 
tem dielectric function 
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where u p i t i = ^Aire 2 n/e QO m*, e, = Jk B Ti/m*, and 
71 is the carrier densities. Equation (|6]) contains both 
phonon and electron-hole plasma contributions; these are 
the second and third terms, respectively. The LO phonon 
contribution is dispersionless. The plasmon contribution 
contains non-zero wave vector modes which account for 
the longitudinal inhomogeneity of the electron-hole gas; 
this term is. equivalent to the Debye-Huckel screening 
(classical) .E^l The zeros of the dielectric function define 
the normal modes of the system: e(q,u>) = =>• 

Figure 2(a) displays the coupled carrier-phonon modes 
(solid lines) of the system in the long wavelength limit 
(q = 0) as a function of the carrier density. The branches 
L + and L_ arise from the anticrossing of the LD and 
plasmon modes (dashed line) due to coupling.O Note 
that for increasing n the L_ branch tends to the screened 
LO phonon frequency i>to- Note that screening is the 
strongest for q = and high densities [Eq. (0)]. In the 
low-density limit, in contrast, L + converges to the bare 
LO frequency islo- 

Figure 2(b) shows the wave- vector dependence of three 
particular L_ frequencies corresponding to n = 2, 4, and 
10 x 10 17 cm -3 . Here as q increases the frequency of 
these three modes tend to the bare vlo value. This in- 
teresting feature is dua.to screening ineffectiveness for 
non-zero wave vectors, 113 i.e., for large q values screening 
by the electron-hole gas becomes ineffective. Equation 
(||) clearly shows the electron- hole contribution to e(q, to) 
gets suppressed for increasing q. In addition, we find that 
a more rigorous description of screening within the Ran- 
dom Phase Approximation corroborates our semiclassical 
normal- mode analysis [see dashed lines in Eig. 2(b)]. We 
consider only the static limit for simplicity!!^ In this limit 
and for the high temperatures used, the RPA calculation 
is straightforward. Ell 

The mode analysis above makes it clear that the 
physics behind the bare LO phonon mode in our numer- 
ical results (particularly at high densities, Fig. 1) is that 
of screening ineffectiveness. We believe this peak at vlo 
is an L_ unscreened q =^ mode; not a low-density L + 
mode with q = 0. Note that this q ^ phonon-like mode 
is not affected by the strong damping of the electron-hole 
gas (r e (^) = 10 fs for Tep) its damping constant is ac- 
tually that iff the long-lived phonon modes (we assume 
T p h = 4 ps.H) Interestingly enough, the "dip" at i>to in 
Fig. 1 follows from e(q, lj) — > oo as lj — > loto which 
leads to a suppression of the internal electric field. This 
suppression is not complete in our simulation because of 
damping. 

Electro- optic modulation signal. The relevant "probe 
signal" in a pump-and-probe electro-optic modulation ex- 
periment is the quantity S(t) oc Jj J a e~ ax [yE(x, t) + 
@w(x, t)]dx, where (3 and 7 are (electro-optic) constants.El 



Figure 3 shows the jEpurier spectrum of S(t) for an n- 
doped GaAs systemliS with an arbitrarily-chosen small 
absorption length a -1 = 10 nm and no depletion layer 
electric field.Ea In this case longitudinal inhomogeneity is 
relevant. Similarly to the THz emission signal, here we 
also clearly see a peak near the bare LO phonon modo. 
This feature is again due to screening ineffectivenessE2l 
[see Fig. 2(b)]. Note that the reflectivity spectrum is 
qualitatively different from the emission one. No dip at 
the screened LO frequency vto should be expected in the 
reflectivity spectrum when both the phonon and internal 
electric fields contribute to the response (1/3/7 1 > 1 hi 
GaAs). In this case, the response does not vanish when 
E(x, t) is suppressed. In our model, no lateral averaging 
over the pump spot is needed in order to obtain peaks at 
vto and v LO - 

Lateral vs longitudinal inhomogeneity. A pump-and- 
probe reflectivity measurement r-with an excitation spot 
much larger than the probe onetil (for Te a cy geometry 
should be used) can conclusively discern contributions 
from these distinctive inhomogeneities. By focusing the 
probe right on the high-density central region of the ex- 
citation spot, one should see both bare and screened LO 
phonon peaks in the probe spectral response. While the 
screened LO peak would represent the L_ mode with 
q = [Fig. 2(a)], the bare LO peak should follow from 
screening ineffectiveness due to q ^ modes [Fig. 2(b)]. 
The presence of this bare phonon mode at high densities 
would constitute a unique signature for the relevance of 
the longitudinal-inhomogeneity effect we propose in this 
paper. Multiple quantum wells and superlattices are also 
suitable systems to search for this effect;c3 the periodic 
modulations along the growth direction give rise to lon- 
gitudinal inhomogeneity. Recent reflectivity experiments 
in GaAs quantum-confined geometriesu have found both 
screened and bare LO modes; these experiments, how- 
ever, use comparable pump and probe spot sizes. 

We have studied the coupled dynamics of electrons, 
holes and phonons within a hydrodynamic approach. In 
longitudinally inhomogeneous systems such as Te, we find 
that screening ineffectiveness at non-zero wave vectors 
allows for oscillations at the bare LO phonon frequency 
in both reflectivity and terahertz emission spectra. We 
believe this mechanism provides an appealingly simple 
(alternate/complementary) explanation for the existence 
of these bare LO modes even at high densities. 
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FIG. 1: Calculated THz signals for differing excitation densi- 
ties, in time (a) and frequency domains (b), (c), and (d). The 
emitted spectra (b) display features at lower frequencies due 
to electron-hole separation and at the characteristic phonon 
frequencies vto and vlo- The theoretical detected spectra 
(c) is obtained by convoluting the emitted signals and the an- 
tenna response [dotted line in (b)]. The calculated (c) and 
experimental (d) detected spectra are very similar: for in- 
creasing excitation densities they both present a strong broad 
band followed by a dip at the TO phonon frequency and a 
weaker peak at the bare LO phonon frequency. Note also 
that the signals get stronger for increasing excitation densi- 
ties. Observe that the relative intensities of the broad peak 
and that of the LO peak are reversed after convolution [cf. (b) 
and (c)]. Surprisingly enough, a bare LO mode survives even 
at high densities n ~ 5 x 10 18 cm -3 ; in our model this follows 
from screening ineffectiveness for non-zero wave vectors aris- 
ing from the intrinsic electron-hole longitudinal inhomogeneity 
in Te. This effect does not qualitatively change with T e /Th 
[see inset in (c)]. The inset in (d) shows the experimental 
time-dependent detected signal. 



FIG. 2: Electron-hole-phonon coupled modes. The q = 
modes (a) show the usual two branches L_ and L+ due to 
anticrossing. Here screening is operative thus making vlo —> 
vto for high densities. Screening ineffectiveness is clearly 
seen in (b) for non-zero wave vector modes: as q increases the 
L- branch approaches the bare vlo frequency. Both Debye- 
Hiickel and RPA show this trend, which is consistent with our 
numerical findings. 



FIG. 3: Calculated reflectivity spectra for an n-doped system 
with short absorption length. Illustratively we use GaAs with 
an arbitrarily small a -1 . Note that both vto and vto are 
present in the response. No averaging over the spot profile is 
performed. 
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